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Abstract 

We calculate perturbatively the effect of density on hadronic prop- 
erties using the chiral quark model implemented by the QCD trace 



anomaly to see the possibility of constructing Lorentz invariant La- 
grangian at finite density. We calculate the density dependent masses 
of the constituent quark, the scalar field and the pion in one-loop 
order using the technique of thermo field dynamics. In the chiral 
limit, the pion remains massless at finite density. It is found that the 
tadpole type corrections lead to the decreasing masses with increas- 
ing baryon density, while the radiative corrections induce Lorentz- 
symmetry-breaking terms. We found in the large Nc limit with large 
scalar mass that the tadpoles dominate and the mean-field approxima- 
tion is reliable, giving rise a Lorentz-invariant Lagrangian with masses 
decreasing as the baryon density increases. 

PACS numbers 21.65.+f, 25.75.+r, 12.39. Fe Ki, 24.85.+p 



1 Introduction 

It has been discussed by a number of authors that QCD trace anomaly 
plays an important role at finite density and temperature. Dilaton fields 
which mimic the QCD trace anomaly are used to scale non-linear cr-model to 
study phase transitions and Nambu-Jona-Lasinio model to study quark 
and gluon condensate [Q] at finite density. A relativistic hadronic model with 
dilaton fields for nuclear matter and finite nuclei is constructed and applied 
to the one-baryon loop level to finite nuclei [Q. It has been also suggested by 
Brown and Rho that the scale anomaly of QCD suitably incorporated into 
an effective chiral Lagrangian allows one to study dense-medium effects in a 
simple mean-field approximation provided the parameters of the Lagrangian 
scaled in medium. The predicted scaling law (called "BR scaling" ) has been 
found to be consistent with a variety of nuclear phenomena ^ • Starting 
with chiral quark model of constituent quarks and pions that incorporates 
the QCD scale anomaly, Beane and van Kolck|p showed that the "mended 
symmetry" of Weinberg can be realized in the dilaton limit with the ef- 
fective Lagrangian. As discussed in 0, there is a close relation between the 
physics of BR scaling and that of the dilated chiral quark model of Beane 
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and van Kolck 0. Thus the dilated chiral quark model is supposed to de- 
scribe in some qualitative and heuristic way the state of matter near the 
chiral phase transition much as BR scaling is to describe the "dropping" 



hadron masses near the phase transition In the recent work[|TT[], this 

aspect of the dilated chiral quark model was exploited to predict the tem- 
perature dependence of hadronic properties near the phase transition: this 
calculation confirms what has always been tacitly assumed, namely, a more 
rapid decreasing of f-j^, m^, and m with temperature in the presence of matter 
density compared with the zero-density a model. In these developments, the 
question as to how relevant the dilaton limit is to dense nuclear matter was 
not, however, properly addressed: for example, how the Lorentz-invariant 
description of dilaton limit in the presence of Fermi sea could be a good 
approximation. 

In this paper, we make an attempt to answer the question by calculating 
density effects starting from the chiral Lagrangian implemented with the 
QCD trace anomaly. We calculate perturbatively the effect of density on 
hadronic properties using the chiral quark model with dilaton fields from 
which Bean and van Kolck derive, in the dilaton limit, the dilated chiral 
quark model. Since we are interested mostly in density effects, we take the 



low-temperature limit and focus on the parts that explicitly depend on the 
density. Using large approximation [Q, we show that the terms which 
break Lorentz invariance can be neglected at large Nc and for large dilaton 
mass and discuss the physical relevance of Lorenz-invariant Lagrangian at 
finite density. 

The paper is structured as follows. In section 2, we present the chiral 
quark model with the dilaton fields and analyze the one-loop structures of 
the model using large approximation. We calculate the density-dependent 
masses of the pion (m^), the constituent quark (m) and the scalar field {m^) 
in section 3 at one-loop order using the technique of thermo field dynam- 
ics. In section 4, we discuss the relevance of mean-field approximation in 
our perturbative calculations and summarize the results. The basics of the 
thermo field dynamics are summarized in Appendix I. Some details of our 
calculations are given in Appendix II. 

2 Effective Chiral Lagrangian 

The effective chiral-quark Lagrangian implemented with the QCD 
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conformal anomaly is^] 

— — TTl — 1 f'^ 

L = f)^ + 9A^ Al^i^ - ^i^i^X + -^tr{d,Ud''U^)x' 

Id ^Jd 

+\d,xd^X-V{x) + - (1) 

where = ^(^^9^^ + id^^), = \i{ed^.i - id^^) with e = U = 
exp(^^j^) and qa = 0.752 is the axial- vector couphng constant for the con- 
stituent quark ^ . The potential term for the dilaton fields, 



nx) = -^bx^-X^ln(f,)], (2) 
^Jd ^ Id 

reproduces the trace anomaly of QCD at the effective Lagrangian level. We 
assume that this potential chooses the "vacuum" of the broken chiral and 
scale symmetries, < 0|x|0 >= fd and the dilaton mass is determined by 

'"x ~ dx^ \x=fd- 

After shifting the field, X ^ fd + x! ■, eq-(0) becomes 



— — — Tfl — 

fd 



^ Since the role of gluons in the chiral quark model is negligible, we will ignore the terms 



containing gluons|13| |14| . 

^The scalar field that figures in the dilaton limit must be the quarkonium component 

that enters in the trace anomaly, not the gluonium component that remains "stiff' against 

the chiral phase transition. See ||] for a discussion on this point. 



+ ]f^trid,Ud^U^) + lktrid,Ud^U^)x' + ]^trid,Ud>^U^)x" 
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+ T^d.x'd'-x' -V{fd + X') + - ■ (3) 



Expanding U in terms of pion fields and collecting the terms relevant for 
one-loop corrections, the Lagrangian for the fluctuating field can be written 



as 



Jd Jn ^/tt 



+ Id.nd'-Tr+^d.nd^nx' + ^d.nd^nx 
^ Id ^Jd 

1. , 1 o ,0 Sm? 



/2 



+ + (4) 

We analyze the Lagrangian using the large A''^ approximation at one-loop 
order. From 't Hooft and Witten[0, we know the Nc dependence of the 



properties of mesons and baryons. Some of them, which are relevant to us, 
are summarized as follows: Three meson vertiex is of order Baryon- 
meson-baryon vertex is of order ^/Wl. But, since the constituent quarks in 
eq.(|l|) have color indices, there should be some modification to the meson- 
fermion vertex. The mesons are created and annihilated by quark bilinears, 

B = Nc^i'"'^" (5) 

If we consider matrix element of eq.(H) between two color-singlet baryons. 



the order is Nc/Ni since N,, quarks are involved to be annihilated by B. 
However in our case the ground-states are the constituent quarks with defi- 
nite color and only one quark which has same color with ground-state quarks 

can contribute the matrix element; thus the quark-meson-quark vetrex is of 
1 

order 1/Nc ■ So quark-meson-quark vertex is of order while there are 
no changes in N,. counting for mesonic vertices. Now we can determine Nc 
dependences of one-loop diagrams. The tadpole contribution to consitutent 
quark mass, Fig. 2, is of order Nc{-^^Y — 1 and that from Fig. 3 is of order 
{^^Y — The two diagrams, Fig. 4 and Fig. 5, for dilaton mass correc- 
tions are of the same Nc order. 

3 Perturbative calculations 

We compute the diagrams representing mass corrections at one-loop ap- 
proximation. At very low temperature, the Bose-Einstein distribution func- 
tion nB{k) — ^3FZT goes to zero but the Fermi-Dirac distribution function 
'npip) = e/^Cp-Mj+i becomes 9{fj.—p). Thus in our case, the finite-temperature 
and -density corrections come from quark propagators. In our calculations of 

7 



Feynman integrals, we shall follow the method of Niemi and Semenoff [|T^ . 
In the case of renormalization of finite-temperature QED one encoun- 
ters a new singularity like ^ ^i3k_i - However in our case, because pion-quark 
vertices depend on momentum, we do not have any additional infrared singu- 
larities at finite temperature and density. At finite temperature and density, 
the lack of explicit Lorentz invariance causes some ambiguities in defining 
renormalized masses. The standard practice is to define density-dependent 
(or T-dependent) mass corrections as the energy of the particle at p = 



1^. In our case, this definition may not be the most suitable one because 
of the Fermi blocking from the fermions inside the Fermi sphere. We find it 
simplest and most convenient to define the mass at p = 0. 



Pion mass 



We first consider pion mass in baryonic matter at low temperature. There 
are three density- dependent diagrams. Fig. 1, that may contribute to pion 
mass corrections. The first one. Fig. la, vanishes identically due to isospin 
symmetry. Explicit calculation of the two diagrams in Fig. lb gives, for 
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/? ^ oo, 

With the change of variable on the second term, p+k — > k, the above integral 
can be rewritten as 



p ■ A;(p^ - 2A; ■ p) + 2mV 2 2^-2/) 

H 5 77] \o{k -m )sin 

P'^ — 2k ■ p 

d^k . 2m2 2m2 , 



J (27r)3 p2 _)_ . j9 p'^-2k-p 
X 6{k^ - m^) sin^ko- (7) 

Since the pion self-energy is proportional to p^, the pole of the pion prop- 
agator does not change. Hence the pion remains massless as long as there 
is no explicit chiral symmetry breaking. Naively we might expect that the 
pion may acquire a dynamical mass due to dynamical screening from thermal 
particles or particle density, which is the case for QED [0 even with gauge 
invariance. The reason for this is the derivative pion coupling to the quark 
field. In eq.(|^), the terms with p ■ k may give rise to terms which are not 



proportional to p^, which lead to the pion mass correction after dk integra- 
tion. However, those terms are proportional to p"^ + 2k ■ p or — 2k ■ p which 
are cancelled by the their denominators and contribute nothing after the dk 
integration. This feature has been explicitly demonstrated in eq.(0). 

Since the pole of the pion propagator does not change, we can take the 
renormalization point at = 0. Then the integrand itself vanishes, so there 
is no wave-function renormalization for the pion in dense medium. 

Quark mass 

At one-loop order, the quark self-energy is given by the diagrams of Fig. 
2 and Fig. 3. The diagram Fig. 2 gives 



where ps is the scalar density obtained from the fermionic loop with thermal 
propagator. 



Id -mi 



(8) 



Ps 
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where Pp = 1^^ ~ ^"^^ With / defined as ps = imi, we get 

- -('jf'^' (10) 

The radiative correction from the pion field in Fig. 3a is 

^'Q{p)-lifm + m)I]. (11) 
Finally, the radiative correction due to the x-field in Fig. 3b is found to be 

where / is defined 

In sum, the self-energy of the quark with four momentum {E,p) can be 
written in the form 

J:q{E,p) = aE^o + b^-p + c- C^f^I (14) 

Id rn^ 

where a, b, and c are 



^ h Jd P 
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where we have used the relation |T6[, J-7 = '^'"^^'^ . The detailed calculations 
of and J ■ p are given in appendix II. 

From eg ■ ([TSD , we can see that the scalar field radiative corrections violate 
the Lorentz symmetry, that is, a ^ —6, unless J° = -^J ■ p. On the other 
hand, the pion radiative corrections in a and b preserve Lorentz covariance. 
Since we have only 0(3) symmetry in the medium for the quark propagation, 
we will adopt the conventional definition of mass as a zero of the inverse 
propagator with zero momentum. The inverse propagator of the quark is 
given by 

G-\E,p) = ^-m-J:{E,p)) (16) 
= E(l-a)7o-(l + 6)7-p-c-m (17) 

The mass is now defined as the energy which satisfies det{G~^) = with 
f = 0, 

{l-a)E = c + m. (18) 

Since a and c are perturbative corrections, eq.(|l^) can be approximated to 
the leading order as 

E = m- ( — y—I + c + am. (19) 

fd m{ 
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If we define the mass as the energy of the particle at finite density, the quark 
mass becomes 

^ ,m 24m^ 

m = m- [ — ) —I 
fd ml 

The first line in eg. (pO]) is just the result of the mean-field approximation, 
in which only the tadpole diagram. Fig. 2, is taken into account, as further 
elaborated on in the section 4. While the dropping of the quark mass with 
density is obvious in the mean-field approximation, it is not clear whether it is 
still true when the radiative corrections. Fig. 3, are included as in the second 
and the third lines in eq. (|20|) . Thus the result obtained at the one-loop order 
does not indicate in an unambiguous way that the quark mass is scaling in 
medium according to BR scaling 0. The specific behavior depends on the 
strength of the coupling constants involved in the theory, gA,m,fT, and fd- 
This does not seem to be the correct physics for BR scaling as evidenced 
in Nature. However, as discussed in section 2, the tadpole contribution to 
consitutent quark mass, Fig.2, is of order 1 and that from Fig.3 is of order 
So we can neglect the Fig.3 for large Nc and obtain in-medium quark 
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mass, 



m* — m — (— ) 




(21) 



In the large Nc limit, therefore, the quark propagator retains Lorentz covari- 
ance and the mass does decrease according to the mean-field approximation. 

Scalair mciss 

Now consider the mass shift of the dilaton field (see Fig. 4 and Fig. 5) . 
The tadpole diagram. Fig. 4, gives 



where the factor 3! comes from the topology of the diagram. This corresponds 
to the mean-field approximation at one-loop order. 

The analytic expression for the contribution from Fig. 5 can be obtained 




fd -m^ 



(22) 



in the large scalar mass approximation, m? « m? 



'X' 




1 
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— m 



1 
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Note that the terms on the second and the third hnes contribute only to 
the energy of the dilaton field and hence break Lorentz invariance. This 
contribution, (P5|), however can be neglected since it is suppressed compaired 
to that of tadpole, eq. (|22D , by Hence the Lorentz invariance is maintained 
approximately in the large scalar mass approximation. Now the renormalized 
dilaton mass becomes 

m*; = ml- 20(^)^1. (24) 

Jd 

The mass of the dilaton drops as density increases. This is consistent with 
the tendency for the scalar to become dilatonic at large density, providing 
an answer to the question posed at the beginning. 

4 Discussion 

So far, we have not taken into account the fact that introducing the ther- 
mal fiuctuation can cause further shifting of the vacuum expectation value of 
the dilaton field. The vacuum expectation value of x'-held in eq.(^) is zero, 
< 0|x'|0 >=0, at zero density. The introduction of the new ground state. 
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I-F >, which are already occupied by constituent quarks with E<Ef{ieTmi 
energy), imphes the necessity of shifting the vacuum. Imposing the condition 
that all tadpole graphs vanish on the physical vacuum, in terms of Weinberg's 
notation [|l^, the condition is 



1—)X'=<F\X'\F> 



+ T = (25) 



dx' 

where P{x') is a polynomial in x' and T is the sum of all tadpole graphs. In 
our case, we are interested in the small change of vacuum expectation value, 
<F\x \F> ^ So we can drop higher powers of the x'-field and retain only 

Jd 



the tadpoles that depend on density at one- loop order. Then the eq.(|25|) 
reads 

- m^Xo - -TPs = (26) 

Id 

where xo =< -^Ix'l-^ >• Then, we get the vacuum expectation value of the 
x'-field. 

^0 = — rrP' (27) 

m^fd 

This is equivalent to Fig. 1 without the external quark line. With the shift 
of the x'-field around Xo, i-e. x' ^ Xo + X in eq.(|D or x /d + Xo + X in 
eq.(|l|), the Lagrangian is modified effectively to 

— — Tn — 7Tl — 

Id Id 
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+\^ifd + Xoftr{d,Ud>^U^) 

J d 

+\d,xd^X - ^(/d + Xo + X) + • • • (28) 

where the elhpsis stands for the interaction terms including pions and x- 
fields. The Lagrangian, eg. (pS]) , shows exphcitly the density dependence of 
fl, rn? and m: 



m 



Id Id 



dx' /J -^^^ ^ 

mJ(l + ^f(l + 3^)^mJ(l + 5^) 



= ^(/, + xo)=m(l + ^). (29) 
Jd Id 

This is the result obtained in the previous section. Hence we can see that the 

mean field approximation in the chiral quark model coupled to dilaton is a 

good approximation in the large limit with massive dilaton, ^i 

and gives 

K (30) 



as predicted by BR scaling. This is somewhat different from the observa- 



tion of Nambu-Freund model which consists of a matter field ip and a 
dilaton field (p. After spontaneous symmetry breaking the matter field and 
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dilaton field acquire masses. There are universal dependences on the vacuum 
expectation values both for the matter field and the dilaton field. One can 
easily see that the universal dependence on the vacuum expectation value 
is no longer valid in our calculations due to the logarithmic potential from 
QCD trace anomaly. 

In summary, we have found that the tadpole type corrections lead to the 
decreasing masses with increasing baryon density, while the radiative correc- 
tions induce Lorentz-symmetry-breaking terms. The pion remains massless 
at finite density in the chiral limit. In the context of large Nc approximation 
with large scalar mass, tadpoles dominate and the mean-field approximation 
is reliable, giving rise to a Lorentz-invariant Lagrangian with masses de- 
creasing as the baryon density increases according to eq. (|30D . This analysis 
with large A''^ approximation gives a clue to construct the Lorentz invari- 
ant lagrangian, i.e, dilated chiral quark model which incoporate the mended 
symmetry of Weinberg into chiral quark model, at finite density. The dilated 
chiral quark model can emerge in dense and hot medium - if it does at all 
- only through nonperturbative processes (e.g., large Nc expansion) starting 
from a chiral quark Lagrangian. 
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Appendix I 



We summarize the basic notation and definition of Tliermo Field Dynam- 



ics(TFD) |0||2^. TFD is a real time operator formalism in quantum field 
theory at finite temperature. The main feature of the TFD is that thermal 
average of operator A is defined as the expectation value with respect to 
the temperature dependent vacuum, |0(/3) >, which is introduced through 
Bogoliubov transformation. 



< A> = Tr{Ae 



)/Tr{e 



< 0{p) I A I 0{P) > 



(31) 



where /? = H is the total Halmiltonian of the system, and /i is the 
chemical potential. The propagator in TFD can be separated into two parts, 
i.e., the usual Feymann part Gp and density-dependent part Gd, 



Gap — GpaP + GoaP- 



(32) 



For fermion with mass m, 



G 



FafS 



\ 
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Goal) = 2niS{p'^ - m^) + m)^,^ 



( 



sin^ 9p^ \ sin 29p^ 



\ 2 



^ sin 29pQ — siv? 6p^ 



with 



cos 9. 



d{po) 



+ 



0{-Po) 



sin ft 



w (l + e-^)V2 ' (i + e^)V2' 
e--/^9{po) e^/^9{-po) 



po 



(l + e-^)V2 (i + e^)V2^ 



where a and (3 are Dirac indices and x = f5{po — /i). 

The propagator of the scalar particle with mass is given by 

A{p)^Af{p) + Ad{p), 

where 



with 



p'^—mg—ie ) 



— 27115 [p^ — ml 



/ 2 1 ^ 

sinh (j)p^ 2sinh20p(, 



|sinh20po sinh^^pQ ^ 



1 



sinh 0, 



(1 -6-12^1)1/2 



po 



(l_e-lyl)i/2' 
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where y = (3po- 



24 



Appendix II 



In this Appendix II, we evaluate various intergrals in the text. 
1. The contribution for the quark self energy from Fig. 3a(pion radiative 
corrections) is given by 



// J i27rrip -ky 
X {-27i)6{k^ -m^){^ + m){^- }^)-i^T'' sin^ Ot^ 
.3 gA^2 f d^k -{p - kf + 2{p - k) ■ k{^- f^,) - m{p - kf 



J (27r)3 {p-ky 



X(5(A;^ -m2)sin2 9, 



ko 



where k denotes \k\. This gives eq.(|ll]) in section 4. 

2. We evaluate J" and J - pin eq. (plSf) , which are needed in the evaluation of 
the quark mass correction from x-field in Fig. 3b. Firstly let us consider J'', 

jO = / -^^^ - m") sin2 Ok, 

J (27r)3 2m2 - 2p ■ - m2 ^ ' 



d?k I 1 .2 

(2^2 2m2-2p-fc-m2 
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^71^ Jo "^^^ 2pk ^""^ -2pk + 2m2 - 2Eko - ml ^ ^^^^ 

Note that due to the factor sin^Okg |/3-*oo^ ^(yU — ko), ko is maximally order 
of chemical potential, i.e., ko ~ /i. Assuming the large dilaton mass, >> 

Q ^ I k 1 1 •2/1 

~ ' (2^2 2m2-2p- A;-m2 '''' 

^ ^ (/x2-,ri2)3/2 (40) 



12tv'^ ml 

A 



Similarly, 



j .p= I P- ^ -5(fc2-m2)sin2^fc„ (41) 

^ ' {2-Kf2m^-2p-k-ml ^ J y j 



can be approximated 



1 ^%(^-m)[/i(/i2-m2)3/2 



87r2 mt 



^ m W/i^ - ln( ^ ^ ~ )] (42) 

3. The contribution from Fig. 5 for the dilaton self-energy, eg. (p3D , is given 

by 

-zSf(/) = -(^)V/^(^+^ + m)(^ + m)(-27r) 
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= (2^[ p^ + 2p.fc + p^-2p.A; ^^^^ -"^^"^ ^'^^ 

„ ,m,r, f dkodx- -r, Eko — pkx + 2171^ 
— 8i{—) / kdk' 



fd J (27r)3 'p^ + 2Eko-2pkx 



,2 



-Eko+pkx + 2m\ 2 2^ ■ 2 n 
^ V^2Ek-^^^^'' 

v"^^2 1 Z'^, rro , P^-4m2 + 2Eko-2pk 
^ ~'^TJ ^ i '^^"^f' + ^"^ p^ + 2^fco + 2p-fe ^ 
- 4m^ ^ y-2Eko + 2pk 
4^^"^ p^-2£;feo-2p-fe ^J 

In the large dilaton mass limit, {-^Y 

Id J P ^P 

TT^ Jd ^ P^ P^ 

= ^2^^ - "^)(/^VA^ - m2 - m /n( ^^^^ ) 

A.m? 8 

-^ln( ^^ + ^^'^'^ )] (44) 
8 m 

where p and k denote \p\ and |^| respectively, and p^ — is used. It is 
suppressed by order of ^ compared to the tadpole contribution. 
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Figure Captions 

Fig.l Self-energy diagrams for pions. Solid lines represent the constituent 
quarks and dashed lines are for pions. D(F) represents thermal(Feynman) 
propagator. 

Fig. 2 Tadpole type diagram for quark. Solid lines represent the constituent 
quarks and wavy lines are dilaton field respectively. 

Fig. 3 Self-energy diagram for quark. 

Fig.4 Tadpole type diagram for dilaton field. 

Fig. 5 Self-energy diagram for dilaton field. 
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